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• ■ Abstract 

O 

• ^ ■ We study geographical effects on the spread of diseases in lattice-embedded scale- 

^^. free networks. The geographical structure is represented by the connecting proba- 

bility of two nodes that is related to the Euclidean distance between them in the 
lattice. By studying the standard Susceptible-Infected model, we found that the 
geographical structure has great influences on the temporal behavior of epidemic 

^ \ outbreaks and the propagation in the underlying network: the more geographically 

t^ I constrained the network is, the more smoothly the epidemic spreads, which is differ- 

ent from the clearly hierarchical dynamics that the infection pervades the networks 

■>!::;j- . in a progressive cascade across smaller-degree classes in Barabasi-Albert scale-free 

^^ ' networks. 

O 

^^ ' Key words: Epidemic spreading; geographical networks; dynamics of social 

^^ . systems 
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•1— I . 

S^ . 1 Introduction 

M'. 

The classical mathematical approach to disease spreading either ignores the 
population structure or treats populations as distributed in a regular medium 
|1J. However, it has been suggested recently that many social, biological, and 
communication systems possess two universal characters, the small-world ef- 
fect [2] and the scale-free property [3], which can be described by complex 
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networks whose nodes represent individuals and links represent the interac- 
tions among them P]. In view of the wide occurrence of complex networks in 
nature, it is important to study the effects of topological structures on the 
dynamics of epidemic spreading. Pioneering works [5 , 6.7j8j have given some 
valuable insights of that: for small- world networks, there are critical thresh- 
olds below which infectious diseases will eventually die out; on the contrary, 
even infections with low spreading rates will prevail over the entire popula- 
tion in scale-free networks, which radically changes many of the conclusions 
drawn in classic epidemic modelling. Furthermore, it was observed that the 
heterogeneity of a population network in which the epidemic spreads may have 
noticeable effects on the evolution of the epidemic as well as the corresponding 
immunization strategies [SllQllTO] . 

In many real systems, however, individuals are often embedded in a Euclidean 
geographical space and the interactions among them usually depend on their 
spatial distances [TT]. Also, it has been proved that the characteristic distance 
plays a crucial role in the phenomena taking place in the system [T2]|13II14|[T5] . 
Thus, it is natural to study complex networks with geographical properties. 
Rozenfeld et al. considered that the spatial distance can affect the connection 
between nodes and proposed a lattice-embedded scale-free network (LESFN) 
model [12] to account for geographical effects. Based on a natural principle 
of minimizing the total length of links in the system, the scale-free networks 
can be embedded in a Euclidean space without additional external exponents. 
Since distributions of individuals in social networks are always dependent on 
their spatial locations, the influence of geographical structures on epidemic 
spreading is of high importance, but up to now it has rarely been studied 
except for p^UTT] . 

In this paper, we study the standard Susceptible-Infected (SI) model on the 
LESFNs, trying to understand how the geographical structure affects the dy- 
namical process of epidemic spreading. Especially, we consider the temporal 
behavior of epidemic outbreaks and the propagation in the underlying net- 
work. It is found that the geographical structure plays an important role in 
the epidemic outbreaks and the propagation of diseases. 



2 The model 



According to Ref. [12], the LESFN is generated as follows: (i) a lattice of size 
N = L X L with periodic boundary conditions is assumed, upon which the 
network will be embedded; (ii) for each site, a preset degree k is assigned 
taken from the scale-free distribution, P{k) ~ k~'^, m < k < K; (iii) a node 
(say i, with degree ki) is picked out randomly and connected to its closest 
neighbors until its degree quota ki is realized or until all sites up to a distance. 
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Fig. 1. (Color online) Degree distribution of an LESFN model with N = 10000, 
7 = 2.5, and for different values of A. The territory parameter A controls the 
influence of the geographical distance on the network structure. 

r(fcj) = A^/k'i, have been explored. Duplicate connections are avoided. Here, 
r{ki) is the spatial distance on a Euclidean plane, denoting the characteristic 
radius of the region that node i can almost freely connect to the others; (iv) 
this process is repeated throughout all the sites on the lattice. Following this 
method, networks with 7 > 2 can be successfully embedded up to a (Eu- 
clidean) distance r(fc), which can be made as large as desired depending upon 
the change of the territory parameter A. Especially, the model turns out to be 
a randomly connected scale-free network when A —^ 00 p^ . Typical networks 
with 7 = 2.5 resulting from the embedding method are illustrated in Fig. [T] 
The power-low degree distributions of the LESFNs achieve their natural cutoff 
lengths for A = 2, 3, and 9, respectively. While for A = 1, the network ends 
at some finite-size cutoff length. 

In order to study the dynamical evolution of epidemic outbreaks, we focus on 
the standard SI model P]. In this model, individuals have only two discrete 
states: susceptible (or healthy) and infected. Each individual is represented by 
a vertex of the network and the links are the connections between individu- 
als along which the infection may spread. There are initially a number of Iq 
infected nodes and any infected node can pass the disease to its susceptible 
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Fig. 2. (Color online) Density of infected individuals versus MC time in the LESFNs 
for different values of A and A. (a) ^ = 9, A = 0.001 (solid line), 0.002 (dash line), 
0.003 (dot line), and 0.005 (dash-dot line), respectively, (b) A = 0.001, A = l (solid 
line), 2 (dash line), 3 (dot line), and 4 (dash-dot line), respectively. The network 
size is A^ = 10000. All the plots were averaged over 500 experiments. 

neighbors at a spreading rate A. Once a susceptible node is infected, it remains 
in this state. The total population (the size of the network) A^ is assumed to 
be constant and, if S{t) and I(t) are the numbers of susceptible and infected 
individuals at time t, respectively, then A^ = S{t) + I{t). In spite of its simplic- 
ity, the SI model is a good approximation for studying the early dynamics of 
epidemic outbreaks [TpO] and assessing the effects of the underlying topologies 
on the spreading dynamics [TSfTg] . 



3 Simulation results 



We have performed Monte-Carlo (MC) simulations of the SI model with 
synchronously updating on the LESFNs. With relevance to empirical evi- 
dence that many networks are characterized by a power-low distribution with 
2 < 7 < 3 |1], we set 7 = 2.5 in the present work. Initially, we select one 
node randomly and assume it is infected. The disease will spread throughout 
the network and the dynamical process is controlled by the topology of the 
underlying network. 

In Fig. [2], we report the temporal behavior of outbreaks in the LESFNs. The 
density of infected individuals, p(t) (= I(t)/N), is computed over 500 real- 
izations of the dynamics. Consistent with the definition of the model, all the 
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Fig. 3. (Color online) Linear-log plots of spreading velocity in the LESFNs for 
different values of A and A. The network size is A^ = 10000. All the curves were 
averaged over 500 realizations. The symbols are the same as in Fig. O 

individuals will be infected in the long-time limit, i.e., limi^oop(^) = 1, since 
infected individuals remain unchanged during the evolution. In Fig. [2]^a), for 
the given connecting region {A = 9), the greater the spreading rate A is, the 
more quickly the infection spreads. Fig. EJ^b) shows the dependence of the in- 
fectious prevalence p{t) on the territory parameter A, while the spreading rate 
is fixed at A = 0.001. As A increases from 1 to 9, the average spatial length 
of edges increases p2]; namely, the nodes have larger probabilities to connect 
to more other nodes, which therefore leads to a faster spread of the infection. 

To better understand the virus propagation in the population, we study in 
detail the spreading velocity, written as [20] 

V,nfit) = ^^p{t)-p{t-l). (1) 



We count the number of newly infected nodes at each time step and report 
the spreading velocity in Fig. [31 Apparently, the spreading velocity goes up to 
a peak quickly and leaves very short response time for us to develop control 
measures. Before the outbreaks of the infection, the number of infected indi- 
viduals is very small, lasting for a very long time during the outbreak, and the 
number of susceptible individuals is very small. Thus, when t is very small (or 
very large), the spreading velocity is close to zero. In the case oi A = 9 (Fig. 
[3]^a)), all plots show an exponential decay in the long-time propagation. These 
results are weakened for a small circular connecting probability (in particular, 
for A = 1), as shown in Fig. [31(b), where the disease spreads in a relatively 
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Fig. 4. (Color online) Temporal behavior of the average degree of the newly infected 
nodes for the SI model outbreaks in the LESFNs with N = 10000 and A = 0.001. 
The data are averaged over 500 dynamical processes. 

low velocity, which therefore slows down the outbreaks. 

A more precise characterization of the epidemic diffusion through such a net- 
work can be achieved by studying the average degree of the newly infected 
individuals at time t [10], given by 



(hnfit)) 



Efcfc4(t)-EfcA:4(t-l) 
lit) - I{t - 1) 



(2) 



where Ikif) is the number of infected nodes with degree k. In Fig. HI we plot the 
temporal behavior of kinf{t) for the SI model with A = 0.001 in the LESFNs. 
Different from the clearly hierarchical dynamics in Barabasi- Albert scale-free 
networks [10], in which the infection pervades the networks in a progressive 
cascade across smaller- degree classes, epidemic spreads slowly from higher- 
to lower-degree classes in the LESFNs. The smaller the territory parameter 
A is, the more smoothly the infection propagates. Especially, strong global 
oscillations arise in the initial time region, which implies that the geographical 
structure plays an important role in early epidemic spreading, independent of 
the degrees of the infected nodes. 
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Fig. 5. (Color online) Temporal behavior of the variance of the density of infected 
nodes in the LESFNs for different values of A. The results are obtained for A = 0.001 
and on networks of size N = 10000. 

Since the intrinsic stochasticity of the epidemic spreading makes each realiza- 
tion unique p31IT7] , it is valuable to analyze the statistical fluctuations around 
the average behavior for assessing simulation results with respect to real out- 
breaks. We measure the variance of the prevalence, defined by 



Var{p{t)) = E[{p{t)-Xt)y]. 



(3) 



In order to evaluate this quantity, we have performed 500 independent runs 
with different configurations of the intrinsic frequencies as well as different net- 
work realizations. Fig. O displays the time evolution of Var{p{t)) for LESFNs 
with different values of A. Since the initial prevalence is fixed and is the same 
for all instances, Var{p{t)) is initially equal to zero and can only increase. At 
very large time, almost all nodes are infected, implying that linif^oo Var{p{t)) = 
0. Compared with Fig. [3I^b), one can easily find that time regions in which the 
fluctuations are maximal are the same as that in which the spreading veloci- 
ties are the fastest. Moreover, an important feature is the rough symmetry of 
curves, regardless oi A = 9 (more randomly scale-free) or A = 1 (more locally 
constrained) and the time regimes in which the fluctuations are maximal cor- 
responding to a small diversity of the degrees of the infected nodes. This is 



also different from that in Barabasi- Albert scale-free networks [TSl. 



4 Conclusions 



We have studied geographical effects on the spreading phenomena in lattice- 
embedded scale-free networks, in which a territory parameter A controls the 
influence of the geographical structure on the network. We have investigated 
the temporal behavior of epidemic outbreaks and found that the spreading 
velocity reaches a peak very quickly in the initial infection period and then 
decays in an approximately exponential form in a more randomly scale-free 
region, which is consistent with previous studies. While the networks are more 
graphically constrained, this feature is apparently weakened. Furthermore, we 
have studied the propagation of the infection through different degree classes 
in the networks. Different from the clearly hierarchical dynamics in which the 
infection pervades the networks in a progressive cascade across smaller- degree 
classes in Barabasi- Albert scale- free networks, epidemic smoothly spreads from 
higher- to lower-degree classes in the LESFNs. Finally, we have analyzed the 
prevalence fluctuations around the average epidemic process. A rough sym- 
metry of curves is found and the time regions in which the fluctuations are 
maximal correspond to a small diversity of the degrees of the infected nodes, 
which is also different from that observed from Barabasi- Albert scale-free net- 
works. 
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